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Abstract 

Near-field  microwave  nondestructive  evaluation  (NDE)  techniques  have  shown  great  potential  for  disbond  de¬ 
tection  in  multi-layer  dielectric  composite  structures.  The  high  detection  capability  associated  with  these  techniques 
stems  from  the  fact  that  near-field  microwave  signals  are  sensitive  to  minute  variations  in  the  dielectric  properties 
and  geometry  of  the  medium  in  which  they  propagate.  In  the  past,  the  sensitivity  of  the  near-field  microwave  NDE 
techniques  to  the  presence  and  properties  of  disbonds  in  multi-layer  dielectric  composites  has  been  investigated  ex¬ 
tensively.  However,  a  quantitative  disbond  thickness  estimation  method  has  yet  to  be  introduced.  In  this  paper,  we 
propose  a  maximum-likelihood  (ML)  disbond  thickness  evaluation  method  utilizing  multiple  independent  measure¬ 
ments  obtained  at  different  frequencies.  We  also  introduce  a  statistical  lower  limit  on  the  thickness  resolution  based 
on  the  mean  squared  error  (MSE)  in  thickness  estimation  and  a  given  confidence  interval.  The  effectiveness  of  the 
proposed  ML  method  is  also  verified  by  comparing  simulation  results  with  actual  measurements. 
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I.  Introduction 

QUANTITATIVE  nondestructive  evaluation  (NDE)  has  significant  utility  in  an  array  of  applications. 

One  particular  application  where  quantitative  NDE  is  frequently  employed  is  in  the  detection  and 
evaluation  of  planner  disbonds/delaminations  in  multi-layer  dielectric  composite  structures.  In  this  context, 
the  NDE  modality  must  be  capable  of  detecting  the  disbond  as  well  as  closely  evaluating  its  thickness.  To 
this  end,  near-field  microwave  NDE  techniques  utilizing  open-ended  rectangular  waveguides  have  shown 
great  promise  [l]-[7]. 

Detection  and  evaluation  of  disbonds  in  multi-layer  dielectric  composite  structures  using  near-field  mi¬ 
crowave  NDE  techniques  have  been  investigated  extensively  in  the  past  (see  [2]  for  example).  The  physical 
model  which  describes  the  interaction  between  microwave  signals  and  a  multi-layer  structure  has  been  de¬ 
veloped  and  validated  experimentally  [2],  [3],  The  sensitivity  of  the  phase  and  magnitude  of  the  reflection 
coefficient,  calculated  or  measured  at  the  aperture  of  the  open-ended  waveguide  probe,  to  variations  in  a 
disbond  thickness  was  demonstrated  in  [5],  In  [4],  a  similar  model-based  approach  was  used  to  estimate 
the  thickness  of  synthetic  rubber  sheets  using  a  root  finding  scheme.  Furthermore,  optimizing  the  measure¬ 
ment  sensitivity  by  selecting  the  appropriate  standoff  distance  and  frequency  of  operation  was  considered 
successfully  in  [5],  [6], 

The  majority  of  the  relevant  work  in  the  literature  emphasizes  the  potential  of  microwave  NDE  techniques 
for  detecting  various  disbonds  and  estimating  their  thicknesses  rather  than  the  estimation  method  or  algo¬ 
rithm  itself.  The  estimation  algorithm  is  mainly  a  functional  procedure  that  works  solely  on  the  measured 
data  to  produce  a  close  estimate  of  the  disbond  thickness.  A  robust  algorithm  should  exhibit  immunity 
against  measurement  uncertainties.  Such  uncertainties  include  system  noise,  intrinsic  errors  in  determining 
the  nominal  values  of  the  dielectric  properties  and  thicknesses  of  the  other  layers  in  the  layered  composite 
structure,  and  the  measurement  system  calibration  errors.  Moreover,  the  estimation  algorithm  should  pro¬ 
duce  unique  (unambiguous)  estimate  of  the  disbond  thickness.  These  requirements  are  crucial  for  robust 
disbond  detection  using  near-field  microwave  NDE  techniques  since  in  general,  the  measurement  parameter, 


3 


i.e.,  the  complex  reflection  coefficient  (T),  is  nonlinearly  related  to  the  disbond  thickness  [2].  The  nonlin¬ 
earity  in  the  measurement  parameters  compounded  with  the  presence  of  measurement  uncertainties  renders 
root  finding  techniques  not  feasible  when  considering  a  relatively  wide  range  of  disbond  thicknesses.  An 
accurate  estimation  algorithm  which  is  capable  of  estimating  unambiguously  the  disbond  thickness  in  the 
presence  measurement  uncertainly,  using  near-field  microwave  techniques  employing  open-ended  rectangu¬ 
lar  waveguides,  has  not  yet  been  introduced  [8]. 

To  the  address  the  above  need,  a  maximum-likelihood  (ML)  disbond  thickness  estimation  algorithm  uti¬ 
lizing  multiple  independent  measurements  obtained  at  different  frequencies  has  been  originally  proposed 
in  [8].  In  this  paper,  we  provide  the  underlying  derivations  for  the  ML  estimator  considered  therein.  Fur¬ 
thermore,  we  extend  the  performance  analysis  and  assessment  of  the  ML  thickness  estimator  beyond  the 
preliminary  investigation  presented  in  [8].  We  also  introduce  a  statistical  lower  limit  on  the  thickness  reso¬ 
lution  based  on  the  mean  squared  error  (MSE)  in  thickness  estimation  and  a  given  confidence  interval.  By 
simulations  and  experiments,  we  show  that  the  proposed  algorithm  produces  highly  accurate  estimate  of  the 
disbond  thickness. 


II.  Maximum-likelihood  Thickness  Estimation 

A.  Measurement  Model 

The  multi-layer  dielectric  composite  structure  configuration  pertaining  to  this  investigation  is  shown  in 
Fig  1.  In  particular,  a  generally  lossy  dielectric  layer  of  relative  permittivity  er  and  thickness  dQ  is  backed 
by  a  conducting  substrate  and  is  irradiated  in  the  near-field  of  an  open-ended  rectangular  waveguide  (RWG). 
An  air-filled  disbond  of  certain  thickness,  cl,  may  be  present  in  between  the  dielectric  layer  and  the  con¬ 
ducting  substrate.  The  objective  is  to  estimate  the  thickness  of  this  disbond,  d,  from  the  measured  reflection 
coefficient  referenced  to  the  waveguide  aperture.  It  is  assumed  that  the  dielectric  properties  as  well  as  the 
thickness  of  the  dielectric  layer  are  known  (i.e.,  er  and  dQ ),  which  is  true  in  a  practical  situation.  We  further 
assume  a  finite  discrete  set  of  disbond  thicknesses  to  be  estimated.  That  is  cl  e  {c/i,  ....  Am]  with  d\  =  0 
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representing  the  no-disbond  case. 

The  open-ended  waveguide  probe  illuminates  the  conductor-backed  structure  with  the  microwave  signal 
and  receives  the  reflected  signal.  At  each  interface  between  any  two  different  materials  in  the  structure, 
the  incident  wave  experiences  a  partial  reflection  due  the  change  in  the  electrical  properties  of  propagation 
medium,  and  undergoes  a  total  reflection  at  the  conducting  substrate.  The  theoretical  model  which  relates 
the  aperture  reflection  coefficient  to  the  structure  constitutive  parameters  can  be  found  in  literature  (see  [2] 
and  [3])  and  it  will  not  be  repeated  here. 

The  complex  reflection  coefficient,  T,  measured  at  the  aperture  of  the  waveguide,  is  used  to  obtain  quanti¬ 
tative  information  about  the  presence  and  thickness  of  the  disbond.  Since  the  same  RWG  is  used  to  transmit 
the  microwave  signal  and  receive  the  reflected  signal  from  the  structure,  the  RWG  radiating  into  the  structure 
can  be  collectively  considered  as  a  single  port  microwave  network  with  the  scattering  parameter  Su  equal 
to  the  aperture  reflection  coefficient,  T.  In  other  words,  the  aperture  reflection  coefficient  can  be  measured 
by  measuring  the  scattering  parameter  Su-  For  this  purpose,  a  calibrated  Shi  type  of  measurement  at  M 
frequencies  is  conducted  using  a  Vector  Network  Analyzer  (VNA) 

Basically,  the  information  sought  about  the  disbond  can  be  acquired  upon  comparing  the  phase  and  mag¬ 
nitude  of  the  measured  complex  reflection  coefficient  to  corresponding  theoretical  values  as  obtained  from 
the  model  mentioned  above.  The  standoff  distance  ds,  the  number  of  frequencies  and  the  frequencies  of 
operation  are  used  as  optimization  parameters  to  maximize  the  estimation  accuracy. 

The  measured  complex  reflection  coefficient  as  a  function  of  frequency,  disbond  thickness,  and  standoff 
distance  can  be  written  as: 

r m(fi,  d ,  ds)  =  ra(fi,  d,  ds)  +  7j,  (1) 

where,  Tm  is  the  measured  complex  reflection  coefficient,  Ta  is  the  actual/theoretical  complex  reflection 
coefficient  as  expected  from  the  model  [2],  ft  is  the  ?’th  frequency  of  operation  fi  £  {/l5  ./*2.  •  ■  •  ,  Im},  d 
is  the  disbond  thickness  to  be  estimated,  ds  is  the  standoff  distance,  and  finally  7i  is  an  uncertainty  term 
that  represents  the  noise  contaminating  the  ith  measurement,  errors  in  the  nominal  values  of  the  dielectric 
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property  and  thickness  of  the  dielectric  layer,  and  the  measurement  system  calibration  errors  .  Henceforth, 
this  term  will  be  collectively  referred  to  as  "noise"  and  it  is  modeled  as  complex  Gaussian  random  variable 
of  zero  mean  and  variance  N0/ 2  per  dimension.  The  parameter  Nq  represents  the  total  power  of  the  noise 
contaminating  the  measurement.  The  choice  of  a  Gaussian  noise  model  is  justified  here  by  the  fact  that  the 
different  sources  of  uncertainties  are  fairly  independent,  and  consequently  their  joint  distribution  function 
tends  to  be  Gaussian  as  suggested  by  the  Central  Limit  Theorem  [9], 

The  objective  of  the  estimation  method  is  to  produce  accurate  estimate  of  the  disbond  thickness  based  on 
the  measured  reflection  coefficient.  Given  the  measurement  model  in  (1),  the  optimum  estimator,  in  Mean 
Squared  Error  (MSE)  sense,  would  be  the  ML  estimator  [10].  The  ML  estimator  is  optimum  as  long  as 
the  disbond  thicknesses  are  equally  probable  in  practice,  i.e.,  the  probability  that  the  disbond  occurs  with 
any  thickness  in  the  finite  set  {c/i,  cZ2,  •  •  • ,  dN}  is  -y.  In  practice,  there  is  rarely  a  priori  knowledge  about 
the  disbond  thickness  being  of  certain  value  with  a  higher  probability  than  any  other  value  in  the  range 
of  interest.  This  situation,  i.e.,  lack  of  knowledge,  is  effectively  modeled  by  considering  all  the  disbond 
thicknesses  as  being  equally  probable.  Otherwise,  the  Bayesian  estimators  outperform  the  ML  estimator 
[10].  Consequently,  for  equally  probable  disbond  thicknesses,  the  ML  estimator  provides  a  lower  limit  on 
the  MSE  performance;  a  limit  against  which  other  estimators  can  be  benchmarked. 

B.  ML  Disbond  Thickness  Estimator 

Since  the  noise  appearing  in  (1)  has  a  Gaussian  distribution,  the  probability  distribution  function  (pdf)  of 
the  measured  reflection  coefficient,  gr,  at  the  / th  frequency  is  given  by: 

1  r  |rm(/i,d,  ds)  -ro(/i,d,  ds)\2 

ttIVo  exp  [  N0 


(2) 


6 


Assuming  that  the  noise  samples,  7^  i  =  1,  2, . . .  M,  are  independent  and  identically  distributed,  i.e.,  being 
Gaussian,  the  joint  distribution  function  of  the  M  measurements  is  found  as: 

M|r,n-ra|2i 

iV0  J  ’  U 

The  log-likelihood  function  when  the  actual  reflection  Ta  is  evaluated  at  the  nth  disbond  thickness,  n  = 
1,  2, . . . ,  N,  is  found  by  taking  the  natural  logarithm  of  both  sides  of  (3)  [1 1],  that  is: 

£(T m;dn)  =  ~M\n(nN0)  -  ^-\Tm  -  Ta |2,  (4) 

iv0 


9r(Tm)  =  7  ,MAfexP 


where  Tm  and  Ta  are  Mxl  vectors  given  by: 


r 


m 


r 


a 


^m(fi  1  d-,  ds^)  d,  C?s) 

r„(/i,d  ni  ds )  r  a(f2,d  ni  ds) 


km(/ M  j  d ,  ds 


£ a[f M 1  dm  ds 


The  maximum-likelihood  estimate  of  the  disbond  thickness,  d  is  the  one  which  maximizes  the  log- 
likelihood  function  given  by  (4).  Since  the  first  term  of  the  right  side  of  (4)  is  independent  of  the  disbond 
thickness,  it  is  sufficient  to  maximize  the  second  term  only.  Based  on  maximizing  the  log-likelihood,  the 
ML  estimator  can  be  formulated  as: 

d  =  min  arg{Z  (dn)},  (5) 

dn 

where  the  ML  decision  metric  Z(di)  is  given  as: 


Z(dn)  = 


(6) 


The  metric  in  (6)  represents  the  square  of  the  Euclidian  distance  between  the  measured  and  actual/simulated 
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reflection  coefficient  vectors.  Basically,  the  ML  estimator  searches  for  the  disbond  thickness  that  minimizes 
this  distance.  This  search  is  conducted  numerically  over  the  disbond  thickness  range  of  interest. 

Since  the  ML  estimator  is  based  on  minimizing  the  Euclidian  distance  between  the  measured  and  ac¬ 
tual/simulated  reflection  coefficients,  estimation  errors  are  most  likely  to  happen  when  the  reflection  coef¬ 
ficients  corresponding  to  different  disbonds  have  close  values  at  the  frequencies  {/i,  /2, . . . ,  /m}-  In  this 
case,  it  becomes  more  difficult  to  differentiate  between  the  different  disbond  thicknesses.  To  reduce  this 
possibility,  the  standoff  distance  and  the  set  of  frequencies  are  selected  such  that  reflection  coefficients  cor¬ 
responding  to  different  disbonds  are  as  distinct  as  possible.  Consequently,  the  standoff  distance  and  the  set 
of  frequencies  that  maximizes  the  estimation  accuracy  of  the  above  estimator  should  be  selected  according 
to  the  following  rule: 

<  F,cls  >=  maxargjy},  (7) 

f,ds 

where,  F  =  {/i,  f2,  ■  ■  ■ ,  /m},  and 

-  M  N 

y  =  nmJ2  E  lr 'a(fi,dn,dj*ra(fi,dk,da)\2,  (8) 

i=  1  n=l,k=l 
k^n 

The  optimization  can  be  accomplished  using  the  theoretical  model  [2]. 

To  inspect  a  structure  similar  to  the  one  described  above  for  disbonds,  the  proposed  estimation  method 
can  be  summarized  as  follows: 

•  decide  on  the  range  of  disbond  thickness  of  interest, 

•  optimize  the  standoff  distance  and  the  set  of  frequencies  to  be  used  according  to  (7), 

•  measure  the  complex  reflection  coefficient  at  the  selected  standoff  distance  and  frequencies,  and 

•  evaluate  the  metric  in  (6)  for  all  possible  disbond  thicknesses. 

Subsequently,  the  disbond  thickness  estimate  would  be  the  one  which  minimizes  (6). 


C.  Statistical  Resolution  Limit 


The  resolution  limit  is  one  of  the  fundamental  performance  metrics  needed  to  assess  the  performance 
of  disbond  estimation  methods.  Basically,  this  limit  establishes  the  minimum  change  in  disbond  thickness 
that  these  methods  can  determine  accurately.  This  is  an  important  practical  issue  which  demonstrates  the 
capabilities  of  thickness  evaluation  methods  in  distinguishing  between  two  close  disbond  thicknesses.  For 
instance,  consider  the  scenario  where  two  structures  are  to  be  inspected.  While  the  first  structure  has  a 
disbond  of  thickness  30  micrometers,  the  disbond  in  the  second  structure  is  45  micrometers  thick.  If  the 
resolution  afforded  by  the  disbond  thickness  estimation  method  is  50  micrometers,  then  this  method  will 
not  be  able  to  distinguish  between  these  two  disbonds.  On  the  other  hand,  the  two  disbonds  can  be  easily 
distinguished  from  one  another  if  the  method  is  capable  of  disbond  thickness  resolution  of  5  micrometers. 
The  resolution  limit  offers  the  needed  practical  insight  into  the  capabilities  of  the  estimation  method  and 
facilitates  performance  optimization,  e.g.,  using  different  frequencies  and  standoff  distances.  It  also  consti¬ 
tutes  a  framework  to  compare  different  disbond  thickness  estimators.  Such  a  resolution  limit  has  not  been 
introduced  previously  for  near-field  microwave -based  inspection  of  disbonds  in  layered  dielectric  composite 
structures.  The  determination  of  the  statistical  resolution  lower  limit  against  which  other  estimators  can  be 
benchmarked  follows. 

The  statistical  nature  of  the  resolution  limit  (as  opposed  to  deterministic)  is  a  direct  consequence  to  the 
random  variations  in  the  measurements  due  to  the  noise.  Therefore,  the  estimated  disbond  thickness  is 
essentially  a  random  variable  with  a  certain  mean  and  standard  deviation.  Since  the  noise  is  zero  mean, 
the  mean  of  the  estimate  is  the  actual  disbond  thickness,  and  its  standard  deviation  is  a  function  of  the 
noise  power.  The  mean  of  estimate  would  wander  around  the  actual  value  of  the  disbond  thickness  in  an 
interval  related  to  the  standard  deviation.  The  confidence  level  that  the  mean  of  the  estimate  is  indeed  in 
the  neighborhood  of  the  actual  disbond  thickness  is  directly  related  to  the  way  we  define  that  neighborhood. 
The  statistical  resolution  limit  we  introduce  herein  is  based  on  defining  that  neighborhood  as  a  function  of 
the  average  mean  squared  error  of  the  estimator  and  the  confidence  level  in  the  estimate. 
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As  mentioned  above,  the  ML  estimator  is  the  optimum  estimator  that  works  on  the  measurement  model 
as  per  (1).  In  other  words,  the  ML  estimator  is  capable  of  providing  the  minimum  MSE  in  estimating  the 
disbond  thickness.  Hence,  the  MSE  performance  offered  by  the  ML  estimator  is  basically  a  lower  limit 
on  the  resolution  for  disbond  thickness  estimation.  The  average  MSE,  over  all  possible  disbonds,  can  be 
expressed  as: 

1  N 

MSE=-J2K-dn\2,  (9) 

71=1 

where  dn  is  the  ML  estimate  of  the  nth  disbond  thickness  dn.  For  p  confidence  interval  of  the  obtained 
estimate,  0  <  p  <  1,  the  resolution  lower  limit,  5,  can  be  expressed  as: 

5  =  2  aVMSE.  (10) 

where  p  is  the  probability  that  the  actual  disbond  thickness  lies  within  resolution  limit  around  the  estimated 
disbond  thickness,  and  the  value  of  a  is  determined  by  solving  the  following  equation. 

1  fX+a  _£2 

p  =  _ _  /  e  2  dx  (11) 

Jx-a 

The  resolution  limit  in  ( 1 0)  is  derived  based  on  the  assumption  that  the  average  estimation  error,  J2n=i  ( dn — 

dn),  is  Gaussian  distributed  with  zero  mean  and  has  a  standard  deviation  equal  to  \/MSE.  This  is  a  reason¬ 
able  assumption  since  the  underlying  noise  is  Gaussian. 

Basically,  the  resolution  limit  given  by  (10)  is  the  confidence  interval  of  the  disbond  thickness  estimate  as 
drawn  from  the  measurements.  Finer  resolutions  can  be  obtained  by  reducing  the  MSE  and/or  the  confidence 
level.  High  confidence  level  indicates  higher  percentage  (from  the  total  number  of  measurement  attempts) 
that  the  actual  disbond  thickness  lies  within  the  resolution  limit  around  the  estimated  disbond  thickness. 
Consequently,  attempting  to  obtain  finer  resolution  by  using  a  low  confidence  level  results  in  lower  percent¬ 
age  that  the  targeted  resolution  is  actually  attained.  The  practical  ramification  of  the  derived  resolution  lower 
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limit  will  be  discussed  in  the  following  section. 


III.  Numerical  Results 

For  simulation,  an  air  disbond  with  a  thickness  varying  from  0  to  0.5  mm,  in  steps  of  0.01  mm,  was 
considered  in  the  X-band  (8.2  -  12.4  GHz)  frequency  range  and  a  standoff  distance  ranging  from  0  to  5  mm 
in  steps  of  1  mm.  The  disbond  was  introduced  under  a  dielectric  slab  with  a  thickness  of  0.778  mm  and  a 
measured  complex  permittivity  of  6.1  —  j'0.37  [12]. 

The  theoretical  complex  reflection  coefficient  was  computed  for  each  combination  of  disbond  thickness 
dn,  standoff  distance  ds,  and  frequency  of  operation  using  the  formulation  mentioned  earlier  [2],  Complex 
Gaussian  noise/uncertainty  with  a  known  power,  N0,  was  added  to  the  computed  reflection  coefficients. 
Thereafter,  the  noisy  reflection  coefficient  was  presented  to  the  proposed  estimation  algorithm  to  estimate 
the  disbond  thickness.  The  signal-to-noise  ratio  (SNR)  in  our  simulation  is  defined  as: 


SNR 


1  yM 

M  i= 1 


r a(fi,d,da)  - 


N0 


(12) 


Starting  with  M  =  21  uniformly  spaced  frequency  points  in  the  X-band  frequency  range,  the  first  step  was 
to  find  the  optimum  standoff  distance  (in  the  set  0, 1, . . . ,  5  mm)  for  detecting  the  disbond  and  accurately 
estimating  its  thickness.  The  (average  over  all  thicknesses)  MSE  as  given  by  (9)  was  used  as  the  figure- 
of-merit  for  this  optimization.  This  is  a  rational  figure  of  merit  since  the  MSE  is  inversely  proportional  to 
selection  metric  in  (8). 

Fig.  2  shows  the  average  MSE  in  estimating  the  disbond  thickness  as  a  function  of  SNR  at  different 
standoff  distances.  It  is  evident  that  the  estimation  accuracy  represented  by  the  MSE  deteriorates  at  low 
SNR,  as  expected.  This  is  attributed  to  the  fact  that  at  low  SNR  the  noise  dominates  over  the  desired  signal, 
i.e.,  the  actual  reflection  coefficient.  As  shown  in  Fig.  2,  the  standoff  distance  of  3  mm  results  in  the  minimum 
MSE  over  the  entire  range  of  SNR  considered.  Hence,  3  mm  is  selected  as  the  optimum  standoff  distance. 

To  illustrate  the  influence  of  the  number  of  frequencies,  M,  on  the  estimator  performance,  the  MSE  with 
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different  number  of  frequencies  is  computed  when  the  standoff  distance  is  set  to  3  mm,  as  depicted  in  Fig. 
3.  It  is  clear  that  the  disbond  estimates  become  less  accurate  as  the  number  of  frequencies  is  reduced.  This 
is  expected  since  by  reducing  the  number  of  frequencies,  less  averaging  is  performed  over  the  measurement 
uncertainty.  Increasing  the  number  of  frequencies  beyond  16,  however,  provides  marginal  performance  gain. 

Based  on  the  results  obtained  above,  the  statistical  resolution  lower  limit  is  computed  as  a  function  confi¬ 
dence  level  for  SNR  of  9  dB  and  18  dB  with  21  frequencies  as  depicted  Fig.  4.  It  is  evident  that  as  the  SNR 
increases,  finer  resolution  can  be  obtained  at  all  confidence  levels.  This  is  mainly  due  to  the  fact  that  the 
MSE  decreases  monotonically  as  a  function  of  SNR  as  observed  in  Fig.  3.  It  is  also  shown  that  the  possible 
resolution  attained  with  higher  confidence  level  is  larger  compared  to  the  lower  confidence  levels.  Basically, 
the  higher  the  confidence  level,  the  wider  the  interval  in  which  the  disbond  thickness  would  exist. 

As  an  example,  consider  the  95%  confidence  interval  (i.e.  p  =  0.95)  curve.  Fig.  4  shows  that  the  minimum 
achievable  resolution  is  around  35  pm  and  13  pm  at  SNR  of  9  dB  and  18  dB,  respectively.  For  SNR  =  9 
dB,  the  resolution  limit  implies  that  if  the  measurements  are  repeated  100  times  and  presented  to  the  ML 
estimator,  95  of  the  times  the  the  actual  disbond  thickness  would  be  in  the  interval  ±17.5  pm  centered 
around  the  produced  estimate  of  the  disbond.  Practically  speaking,  the  ML  estimator,  at  SNR  of  9  dB  and 
95%  confidence  level,  will  not  be  sensitive  to  a  change  of  less  than  35  pm  in  the  disbond  thickness.  If  lower 
resolution  than  35  pm  is  needed,  one  has  to  increase  the  SNR  in  the  system.  For  example,  doubling  the 
SNR  to  18  dB  yields  a  resolution  of  around  13  pm  with  95%  confidence  level.  It  should  be  noted  that,  in 
the  X-band  where  the  maximum  available  bandwidth  is  4.2  GHz,  the  far-held  thickness  resolution  limit  is 
around  3.6  cm  in  air1.  Thus,  the  improvement  obtained  by  using  a  near-held  approach  in  conjunction  with 
ML  algorithm  has  improved  the  attainable  resolution  by  a  factor  of  1200  for  this  example. 

IV.  Experimental  Results 

To  demonstrate  the  efficacy  of  the  proposed  approach,  a  multi-layer  structure  similar  to  the  one  depicted 
in  Fig.  1  was  assembled.  We  used  synthetic  rubber  sheet  of  thickness  4.42  mm  as  the  dielectric  layer.  The 

xThe  far-field  resolution  limit  with  inspection  bandwidth  of  B  is  c/2B  where  c  is  speed  of  light  in  air. 
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dielectric  properties  of  the  rubber  sheet  were  measured  at  X-band  using  two-port  loaded  transmission  line 
technique  [12]  to  be  Er  =  7.2  —  j'0.34.  An  air-filled  disbond  of  varying  thickness  was  introduced  in  between 
the  rubber  sheet  and  the  conducting  substrate  by  moving  the  substrate  away  from  the  rubber  sheet  using 
a  high  precision  positioning  apparatus.  The  disbond  range  of  interest  was  set  to  be  from  0  to  1  mm.  The 
thickness  of  the  disbond  was  varied  in  this  range  with  a  step  of  0.05  mm.  Finally,  the  standoff  distance  was 
fixed  at  3  mm,  and  an  X-band  (WR-90)  open-ended  waveguide  probe  was  used  to  irradiate  the  sample.  In 
practice,  the  disbond  thickness  is  typically  much  smaller  than  the  thickness  of  the  dielectric  layer.  Some  of 
the  values  used  in  this  experiment,  however,  were  for  illustration  purposes. 

At  each  disbond  thickness,  calibrated  swept  frequency  Su  measurements  were  conducted  using  an  HP85 10C 
vector  network  analyzer.  To  randomize  the  measurement  error  (especially  positioning  errors),  the  measure¬ 
ments  where  repeated  6  times  and  for  each  case  we  repositioned  the  substrate  and  the  rubber  sheet.  There¬ 
after,  the  measured  6-sample  data  was  presented  to  the  proposed  algorithm  to  estimate  the  disbond  thickness 
and  the  results  were  averaged. 

Fig.  5  shows  the  computed  metric  from  (6)  as  a  function  of  disbond  thickness  with  different  number  of 
frequencies  for  a  disbond  thickness  of  0.5  mm.  At  each  number  of  frequencies,  the  ML  estimator  selects  a 
disbond  which  corresponds  to  the  minimum  point  on  the  curve  as  indicated  in  (5).  Table  I  summarizes  the 
estimation  results  for  this  example.  As  shown  in  the  table,  as  the  number  of  frequencies  increases,  better 
estimates  of  the  disbond  thickness  are  obtained.  Increasing  the  number  of  frequencies  beyond  15  however, 
did  not  measurably  improve  the  performance  of  the  algorithm. 

To  further  assess  the  capabilities  of  the  proposed  algorithm  for  providing  high  resolution  disbond  thickness 
estimates,  the  estimates  were  computed  for  all  disbonds.  Fig.  6  shows  the  ML  estimates  as  a  function  of  the 
actual  disbond  thickness  with  different  number  of  frequencies.  A  perfect  estimation  trace  is  also  depicted  on 
the  graph  as  a  performance  reference.  It  is  evident  that  the  proposed  algorithm  was  able  to  distinguish  each 
disbond  thickness  independently  and  provide  good  estimate  of  its  thickness  especially  when  large  number  of 
frequencies  were  used.  However,  the  algorithm  produced  relatively  large  estimation  errors  when  the  disbond 
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thickness  was  very  small.  This  is  mainly  attributed  to  the  positioning  errors.  Setting  the  disbond  thickness 
to  zero  (i.e.,  the  reference)  was  problematic  due  to  the  elastic  nature  of  the  rubber.  Although  repositioning 
for  each  sample  was  partially  successful  in  randomizing  the  positioning  errors  for  large  disbond  thicknesses 
(at  the  expense  of  increasing  the  standard  deviation),  it  resulted  in  a  consistent  thickness  bias  for  the  zero 
disbond  thickness  case.  Consequently,  the  ML  estimator  was  sensitive  to  50  //m  change  between  the  large 
disbond  thicknesses  (it  provided  very  close  estimates  for  these  thicknesses),  and  it  failed  to  distinguish  the 
same  50  /jm  change  between  the  zero  and  the  50  //m  disbond  thicknesses. 

Other  sources  of  inevitable  measurement  errors  including  the  uncertainty  in  the  dielectric  properties  of 
the  rubber  sheet  and  the  mismatch  between  the  theoretical  model  and  measurement  (e.g.,  the  theoretical 
model  neglects  the  higher-modes  effect  and  assumes  infinite  waveguide  flange),  are  deemed  to  impact  the 
performance  of  the  ML  estimator  as  well. 

V.  Concluding  Remarks 

A  quantitative  high-resolution  disbond  thickness  method  has  been  presented  in  this  paper.  Based  on 
Maximum  Likelihood  (ML)  approach,  the  proposed  method  utilizes  multiple  independent  measurements 
obtained  at  different  frequencies  to  estimate  the  disbond  thickness.  The  proposed  method  lends  itself  to 
optimization  with  three  degrees  of  freedom,  namely;  the  standoff  distance,  the  number  of  frequencies,  and 
the  frequencies  of  operation. 

By  investigating  the  performance  of  the  ML  estimator,  we  introduced  a  statistical  lower  limit  on  the 
thickness  resolution  offered  by  the  proposed  method.  Since  the  ML  estimator  is  optimum  in  the  MSE 
sense,  the  resolution  lower  limit  provides  a  framework  to  compare  and  analyze  the  performances  of  possibly 
different  estimators. 

By  simulation  and  experiments,  it  was  shown  that  the  proposed  estimation  method  provides  a  promising 
performance  in  detecting  and  evaluating  the  thickness  of  disbonds  in  multi-layer  dielectric  composite  struc¬ 
tures.  We  remark  that  the  proposed  method  might  be  computationally  prohibited  for  real-time  applications 
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where  the  search  for  the  disbond  thickness  should  be  conducted  over  wide  span.  In  such  applications,  it  is 
recommended  to  use  efficient  approximate  ML  estimators. 
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TABLE  I 

Estimation  Results  for  0.5  mm  disbond 


M 

Estimated  Thickness  (mm) 

Percentage  Error  (%) 

6 

0.375 

25 

10 

0.421 

15.84 

15 

0.475 

5 

21 

0.475 

5 

15 
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Fig.  1.  Multi-layer  structure  configuration  with  air- filled  disbond. 
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Fig.  2.  The  MSE  in  estimating  the  thickness  of  the  disbond  as  a  function  of  SNR  at  different  standoff  distances  in  the  X-band  using  21 
frequencies. 


Fig.  5. 


The  computed  metric  from  (  6)  as  function  of  disbond  thickness  with  different  number  of  frequencies. 


Fig.  6.  Estimation  results  for  all  disbond  thicknesses  with  different  number  of  frequencies. 


